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ABSTRACT 



A computer model was developed to simulate a 
projectile ram- propellant feed system for a Liguid 
Propellant Gun. Using a lumped parameter approach, a 
set of si nultaneous differential eguations was derived 
for the complex interaction of the propellant fluid, 
the drivirg injector and the projectile. The computer 
model was verified against a 20 mm experimental 
apparatus. Injector displacement, projectile 
displacement, and chamber pressure were compared for a 
nominal driving pressure of 140 psi. ihe important 
system parameters affecting projectile ram time and 
chamber pressure oscillations were investigated and 
potential prctlem areas for testing with actual 
propellant were identified. 
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I. I N1R0D UCTI QN 



As the state cf the art of weapons technology improves, 
mere complex systems are continually being developed. To 
meet future needs these systems must have improved 
performance, less weight and volume, less complicated 
logistics, and at the same time be more cost effective than 
current systems. One candidate for these future weapon 
systems is the liquid propellant gun (LPG) . Numerous 
invest igatiers over the past thirty years have demonstrated 
its potential if the technology can be developed to utilize 
a liguid propellant in a large caliber gun. 

One impertant area that needs investigation is the 
aspect of fluid propellant handling. The ordnance designer 
must knew the necessary propellant flow rates and pressure 
regimes needed to achieve the desired system performance. 
These parameters must be within allowable safety margins for 
the handling and use of the explosive propellant. One 
possible design of an LPG which would allow a very high 
firing rate utilizes the propellant to ram the loaded 
projectile. Combining the normally separate propellant load 
and ram cycles decreases the required movement of mechanical 
parts. Thus, using this method, it is feasible that firing 
rates in excess cf fpur to five times present rates could be 
achieved. This thesis was directed toward identifying the 
impertant fluid dynamic parameters involved with such a 
projectile ram feed system. 
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II. background 



During 1S74 and early 1975, an investigation was 
conducted at the Naval Postgraduate School to study, both 
analytically and experimentally, the fluid dynamics of a 
liguid propellant under conditions similar to those which 
would exist in a rapid-fire LPG feed system. The results of 
this investigation were to be used to establish such LPG 
design and performance parameters as time-to-load, injection 
supply pressure, injection system configuration, ullage, 
charge-tc-mass ratio, caliber size, and projectile mass. It 
was hoped that this investigation would identify any 
potential problem areas for further detailed research. 

The basic objective of the experimental portion of the 
investigation was to identify what fluid dynamic 
characteristics of a liquid propellant feed system would 
limit leading times and hence rates of fire. A 20 mm 
experimental model of a basic propellant feed system was 
designed and built.' Data on injector displacement , breech 
chamber pressure, and ram gas pressure (input driving 
pressure) , were recorded for driving pressures between 50 
and 220 psig in 10 psi increments. As reported in Ref. 1, 
it was found that the instantaneous behavior of the chamber 
pressure was the result of a complex interaction of inertia 
forces, viscous forces, and the unsteady motion of the 
fluid. The experiments demonstrated that frictional and 
inertial effects were significant during the movements of 
the injector and the projectile slug. Once the projectile 
slug stopped, the effect of entrapped gas in the fluid 
caused large breech chamber pressure oscillations. On 
several runs, sub-atmospheric pressures were experienced 
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which suggested the possibility of cavitation and hence 
vapor-phase ignition-. Over all, it was found that the ram 
tine had a guadratic dependence upon ram pressure. 

The analytical portion of the investigation was directed 
toward predicting the pressure and flew rate of the 
propellant and the projectile slug motion during an LPG 
loading cycle. The mathematical model which was developed 
and programmed on an analog computer predicted the position, 
velocity, and acceleration of the projectile and the 
pressure at various points in the system as functions of 
time. The model was tested against experimental results and 
found to he adequate for the prediciton of projectile ram 
time. An analysis was performed using the model to indicate 
the areas of system redesign likely to be most profitable 
and to obtain preliminary predictions of LEG loading system 
performance under a variety of design conditions. The 
results of the analytical study, as well as a summary of the 
experimental work is contained in Bef. 2. 

The aforementioned model was the starting point for the 
present study. To better understand the derivations in the 
following sections. Figure (1) indicates the geometry of a 
basic projectile ram propellant feed system. An injector 
chamber is filled with propellant, which is then pumped into 
the gun breech by applying ram pressure to the ram side of 
the injector piston. This can be accomplished by using high 
pressure cas from an accumulator or by using an hydraulic 
drive system. The force exerted by the injector piston 
drives the fluid propellant through the connecting line into 
the gun breech. The rapidly accelerated propellant drives 
the already leaded projectile to its seated position ready 
for firirg. 
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Figure 1 - BASIC GEOMETBY 
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The complex nature of this dynamic ram process is not 
easy to describe with mathematical eguations. Many 
researchers have approached the problem using wave mechanics 
and partial differential eguations. These eguations were 
then modified according to the propagation characteristics 
of the system and the various boundary conditions 
encountered therein. This usually led to involved finite 
difference methods of soluticn on a digital computer. 

Another possible approach considers the kinetic energy 
involved in the ram process. This is the approach that was 
used to derive the governing eguations for the previously 
described analog computer model. Beginning with the input 
side, a force balance on the injector piston yielded an 
eguaticn fcr the pressure in the injector chamber in terms 
of the injector motion. Writing Bernoulli's eguations for 
head losses at the expansion to the breech chamber, the 
ccntracticn to the connecting line, and an orifice in the 
connecting line resulted in an eguation fcr the pressure 
drop in the connecting line in terms of the sguare of the 
fluid velocity and tie fluid acceleration. A force balance 
on the projectile slug yielded a third eguation which 
described the breech chamber pressure in terms of the slug 
motion. Tc solve these eguations, seme method of relating 
the motions in the injector, connecting line, and breech 
chamber had to be established. It was assumed that the 
fluid was incompressible and therefore these motions were 
identical. 1c facilitate the analog simulation it was also 
assumed that the input ram pressure was a step input. In 
the process of converting these eguations to an integral 
form for hiring on an analog computer, all variables were 
normalized tc permit scaling. 

Despite these limiting assumptions, the analog computer 
model was able to predict the time to ram the projectile and 
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the chamber pressures which cccur during the time that the 
injector is in motion. Due to the incompr essitility 
assumption, the model was unable to account for the pressure 
transients cr "water hammer" type pressure oscillations 
which occurred when the projectile was seated and the 
propellant fluid was being decelerated. These pressure 
oscillations are important to the system designer because 
the peak breech chamber load pressure is experienced during 
these oscillations. It is also possible that these 
oscillations could interfere with the uniformity cf the 
ignition and subseguent combustion of the propellant. 

h straight forward attempt was made to extend this 
analog model by first converting it to a digital computer 
pregram. This allowed an accurate representation cf the 
input ram pressure and the use of unnormalized variables to 
be incorporated. These improvements increased the accuracy 
of the mcdel but the effect of the compressibility of the 
propellant fluid still was not taken into account. Trying 
to add compressibility effects to this model by adding time 
derivatives cf the pressure terms to the governing eguations 
unlinked the motions of the injector, the connecting line 
fluid, and the projectile slug causing a problem with too 
many unknown variables for the number of eguations involved. 

It was decided that a new approach should be tried which 
would adeguately describe the system pressure oscillations. 
The approach would feature an engineering mcdel which would 
be easy to use, adaptable to any LPG feed system, and not 
obscure the interaction of system parameters by complicated 
mathematics. 
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III. D ERIVAT ION CF COMPUTER MODEL 



A. LUMPEI PARAMETER APPROACH 

The approach that showed the greatest promise for 
modeling the LPG feed system was that of the fluid 
transmission line cpncept. This approach, which has tecome 
popular in the last five to ten years, is based on a 
pressure-voltage and flow velocity-current analogy with 
Electrical Engineering determinations. It is an outgrowth 
of the large amount of effort that has been devoted to 
investigating fluid line transients. Reference 3 is a good 
survey of this field. 

For complex systems this approach is usually simplified 
by using the approximation of lumped parameters. The 
effects of fluid inertia, capacitance, and resistance are 
"lumped" and considered to act only in discreet areas cf the 
system. This results in a reduction cf the unknown 

parameters due tc the assumed lack of interaction cf the 
different fluid effects. This approximation, however, 
results in the necessity of using several empirical 
constants which must be determined by fitting experimental 
data. The ordinary differential equations that are derived 
frcm this method can be solved either by Laplace 

transfer nation or ty computer integration. . 

To consider the effect of fluid inertia it is assumed 
that only pressure and inertia forces are present and that 
• compressibility effects are negligible in the volume under 
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consideraticr . Then, 



dv 

P , - P 2 ■ /° L 7f 



where p is the propellant fluid density, 1 is a 

characteristic length and V =V =V because the flew is 

1 2 

considered irccmpressible for this building block (Fig.2A). 



To consider the effect of fluid capacitance, it is 
assumed that only compressibility effects are important, and 
that inertia and resistance effects may be neglected in the 
volume under consideration. Therefore, 

L dP 

V - V 

1 2 MI 

where K is the effective bulk modulus cf the propellant 
fluid, 1 is a characteristic length and P^=P^=P (Fig. 2B) . 

Eecause fluid resistance can be affected by sc many 
different parameters, it is impossible to write a general 
eguation describing the pressure drop due to fluid 
resistance. It is best to treat it empirically using an 
experimentally derived figure. Thus, in the volume under 
consideraticr , 

P - P = Rv V 

1 2 

where Ev is a function of fluid velocity and V^=V^=V (*ig« 
2C) . Changes in cross sectional area can be accounted for 
by using appropriate area ratios. 

These results specify three building blccks which can be 
combined in any sequence to model the dynamic 
characteristics of a system. The complexity of the model 
can be increased to any degree necessary by including more 
and more cc m hi nations of these three basic building blocks. 
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A) Inertia 
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figure 2 - LOBfED EAE4BETEB BODULES 
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E. LEG SYSTEM GOVERNING EQUATIONS 



Several different combinations of buildirg blocks were 
tried in attempting to model the experimental LPG system. 
The experimental results showed that inertial and resistance 
blocks should be included in the model. The fact that 
pressure transients pccurred in the breech chamber indicated 
that a capacitance block should also be included. To 
minimize the use of computer time, it was decided to start 
with the most simple model thought adequate and work toward 
more complex models as the quality of the computer results 
dictated . 

The first attempt considered that ccmpressibilt y effects 
would dominate in the injector chamber as the accelerating 
piston interacted with the propellant fluid. The resistance 
effects cf changes in the feedline cross sectional diameter 
and associated orifices were modeled as one resistance block 
at the beginning of the connecting line. Inertial effects 
were thought to dominate in the connecting line. This would 
be particularly true in large scale models where the length 
of the connecting line would be very large. It was thought 
that compressibility effects should again dominate in the 
breech chamber as the accelerating propellant fluid drove 
the projectile dcwn the breech chamber. Results in this 
case were governed by inertial effects because the 
calculated connecting line velocity was too great. This led 
to a decrease in the chamber pressure as the projectile 
accelerated until a large negative chamber pressure existed. 
In additicn to this deficiency, the computed system pressure 
cscillaticns were not sufficiently damped. 

To try tc adjust the connecting line velocity, another 
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resistance block was added. As will be seen in the section 
discussing the parameters which affect the pressure 
oscillat icns , the value of this resistance coefficient is 
prcporticral to the amount of damping in the pressure 
oscillations . Even with a value which gave the proper 
damping, the breech chamber pressure was dominated fcy the 
pressure loss due to the accelerating fluid in the 
connecting line. 

To reduce the dominance of the fluid inertia effects, it 
was decided tc restructure the model. To model the injector 
chamber it was decided to consider the effect of fluid 
capacitance anc inertia. The connecting line was modeled by 
twc resistance blocks separated by a fluid capacitance 
blcck. The breech chamber was modeled using a fluid 
capacitance and a fluid inertia block, as was done for the 
injectcr chamber. This model gave the final results 
described in the section comparing computer and experimental 
data. Ey combining these building blocks, the governing 
eguaticns became: (See Fig. (1) for notation) 



For the in“‘ectcr: 



i. Pr A -PA = Mp V + Kfs V + Pdi A 
111 i l i 



2. V - V = P n 

l 2 K i 



3. P - P = cl V 

2 3/12 



For the connecting line: 
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4. P - P = Rv V 

3 4 14 



L2 * 

5. v - V = — p 

4 5 K 4 



6- P - P = Rv V 

5 6 2 5 



For tt€ hreech chamber: 



7. 




6 



8 . 




^3 

K 



P 

8 



9. 



P A - Pds A = Ms V + Kfs V 

8 3 3 8 8 



These equations were subsequently modified tc account 
for changes in cross-sectional area and changing geometry as 
the injectcr chamber decreased in volume and the hreech 
chamber increased. Once the projectile reached the end of 
the breech chamber, the force balance eguaticn was nc longer 
considered applicable. At this time, it was noted that the 
remaining equations for the connecting line and the hreech 
chamber cculd be combined to form a second order 
differential equation for the breech chamber pressure which 
described the "water hammer" type pressure oscillations 
displayed by the experimental data. Thus, 
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K 



K 



" A 3 Rv 2 • K 

10 . P + — — — P + P 

8 A 2 />L 3 8 yOL 3 L 3 8 



p \- 3*- 3 



P 

4 



Once the injector piston reached the end of its travel, 
the injector chamber no longer existed and the governing 
eguaticns were considered as no lcnger applicable. 
Ultimately, the injector chamber pressure was propagated 
down the connecting line as the propellant fluid came to 
rest. Throughout these analyses, input ram pressure was 
modeled using expcrential terms to fit experimental data. 



The eguaticns from the lumped parameter approach fcr the 
propellant fluid were combined with the force balance 
eguations fcr the injector piston and the projectile using 
State Variable methods. Subsequently , a computer program 
was constructed to solve these state variable relations 
using a fourth crder Runga-Kutta integration routine for 
siir ultan ecus first order differential eguations which was 
developed at the Naval Postgraduate School. A brief 
description cf the state variable method and a listing of 
the computer program can be found in Appendix A. 

The fcllcwing is a listing cf the nomenclature and the 
values cf constants used in the derivation cf the governing 
eguations anc the computer program: 



A 1 


Injector 

Cross-Sectional Area 


1.77 


2 

in 


11. 42 


2 

cm 


A2 


Connecting Line 
Cross-Sectional Area 


.255 


2 

in 


1.65 


2 

cm 


A3 


Ereech Chamber 


.49 


2 

in 


3.16 


2 

CE 




Cross-Sectional Area 










A B 


Ram listen 

Cross-Sectional Area 


1.77 


2 

in 


11.42 


2 

CE 


11 


Injector Length 


1.60 


in 


4.22 


CE 
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L2 


Connecting Line Length 


16.0 


in 


40.6 


cm 


L3 


Breech Chamber Length 


5.0 


in 


12.7 


cm 


MF 


Injector Piston Mass 


.0052 


2 

lb-sec 


913. 8 


gm 








in 






MS 


Projectile Mass 


.00053 


2 

lb-sec 


93. 1 


gm 








in 






KFP 


Injector Frictipn 


.001 


lb-sec 


175.7 


gm/sec 




Factor 




in 






KFS 


Projectile Friction 


.001 


lb-sec 


175.7 


gm/sec 




Factor 




in 






FDI 


Injector Back Pressure 


24.0 


PSI 


163.3 


kFa 


FES 


Projectile Back 


.50 


PSI 


3.4 


kFa 




Pressure 










K 


Effective Bulk Modulus 


3200 


PSI 


6808 


kFa 


BBC 


Propellant Density 


.000093 


2 

lb-sec^ 


.994 


3 

gm/cm 








4 

in 






PB 


Bam Pressure 




Computed 




PI 


Injector Chamber Pressure 


Computed 




VI 


Injector Pistcn Velocity 




Computed 




V2 


Injector Chamber Exit Velocity 


Computed 




P4 


Connecting Line Pressure 




Computed 




V6 


Breech Chamber Entrance 


Velocity 


Computed 




F8 


Ereech Chamber pressure 




Computed 




V8 


Projectile Velocity 




Computed 





The treatment cf several areas of the computer model 
continually reappeared as requiring refinement. The 
governing differential equations worked well during the 
dynamic pcrticn cf the feed cycle but experienced 
difficulties during the initial and final static periods. 
The handling cf static friction and back pressure as a 
constant value created the possibility of negative 
velocities urtil the driving ram pressure overcame the 
system tack pressure. These negative values never occurred 
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in the real system due to the geometric restraints or the 
injector piston and projectile; therefore, the computer 
model had to be manipulated to maintain this condition. 
Unfortunately the need for simultaneous solution of the 
governing eguations made it difficult to manipulate the 
initial conditions without greatly increasing the complexity 
of the computer model. After the projectile stopped, the 
transition from the force balance equation on the projectile 
to an equation describing the fluid velocity in response to 
the pressure transients was awkward. No simple differential 
eguation describes the complexity of the wave mechanics 
involved with the reflection of the pressure waves in the 
system. In this regard several alternate approaches 
involving sequential alterations of the governing eguations 
were tried with varying degrees of success. 
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m. description cf exper ime ntal apparatus 



In crdei that the comparison of the experimental and 
computer generated data can be fully understood, a brief 
description cf the NPS experimental apparatus and the 
conduct cf the associated experiments is included (See Ref. 
1 , pages 17-27) . 

The test chamber was fabricated from a three-inch O.D. 
Lucite cjlirder, 18 inches long, bored to a 20 mm inside 
ciameter and fitted with aluminum end caps. The chamber was 
loaded with a trass slug weighing 93 grams which rode cn two 
graphite filled Teflpn sealing rings. The trass slug, which 
simulated the projectile, was cycled from the breech end of 
the chamber tc the barrel end and returned to the treech 
end, completing one hypothetical firing cycle. 

Because cf the desire to vary the charge to mass ratio, 
a variable chamber velocity was necessary. To accomplish 
this with one chamber, a volume control retaining rod was 
designed into the system. This brass rod, bored to allow 
gas to pass its length, was threaded through a plate which 
was attached to the barrel end cap holding the red in the 
chamber. The rod, which has a Teflon disc cn the end, not 
only established the volume of the test chamber, but 
provided a buffer stop for the slug at the end of its 
forward motion. Another Teflon buffer was affixed to the 
breech end cap to cushion the slug in return motion. 

The LEG simulator started a simulated firing cycle with 
the slug at the breech end of the empty chamber, as shewn in 
Fig. (3). 
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Power Piston Injector P —j Test Chamber 
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Figure 3 - PNEUMATIC CONTROL SYSTEM IN RAM POSITION 



The simulated 



This is the ready-to-ram position, 

propellant, distilled water, was then introduced, ramming 
the slug tc the opposite (barrel) end of the chamber as the 
chamber was filled. This was accomplished by applying gas 
pressure tc a power piston which drove the injector piston. 
The injector piston forced the propellant past a flow check 
valve and into the chamber. This placed the slug in the 

ready-tc-f ire cr ir-battery position. In an actual gun, the 
prcpellant wculd he ignited at this time in the cycle. Due 
to laboratory constraints, an expulsion system was used. 

An I VDT displacement transducer was manufactured and 
mounted next tc the injector piston. The LVDT was attached 
tc the connecting rod, between the power piston and the 
injectcr pistcn. The volume of the liquid being placed 
under pressure during each shot was measured by filling the 
system in the ready-to-ram position and then draining it 
into a graduated beaker. By measuring the displacement of 
the injectcr piston head with the LVDT, the volumetric rate 
of fluid injection during the ram stroke was obtained. 

Two pressure taps were drilled in the test chamber. 
These taps were located as close to the breech end as 
possible, cne at 20' degrees frcm top center and the other at 
20 degrees from bottom center in a counterclockwise 
direction, as seen from the breech end of the chamber. Only 
the bottcm lccation was used to record data. 

A 4-chanrel Hewlett-Packard 3960 Magnetic Tape Recorder 
was used to IM record desired data during system operation. 
For each cata run signals were recorded from two Kaman 
diaphragms type ,(1000 psi) pressure transducers, one 

connected to the breech pressure tap on the test chamber, 
and the ether to the gas injection side of the power piston. 
The pressure signals were processed with a Kaman Digi-Vit 
Readout Onit which also provided a visual (digital) display. 
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Data from the IVDT displacement transducer was also recorded 
cn the tape, 

A Brush Recorder (Mark 280) was used to obtain a visual 
display cf the recorded data. By transcribing the desired 
signals cn the Magnetic Tape Becorder at a tape speed cf 15 
feet per minute and playing them back into the Brush 
Becorder at 5-3/4 feet per minute, the time scale of the 
output was expanded by a factor of four on the Brush 
recordings (viz., from a real-time maximum cf 200 mm/sec to 
a delayed tiae maximum of 800 mm/sec.) 
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V 



CCKPARISON OF COMPUTER AND EXPERIM ENTA L DAI A 



As mentioned previously, experimental data was taken for 
ram gas pressures ranging from 50 to 220 psig. Froii this 
collection of data, pne run at 140 psig was selected as 
representative of system performance. It was felt that any 
model which wculd adequately describe systeir operation at 
this intermediate pressure would be valid for the entire 
range of expected IPG driving pressures. 

Figure (4) is a comparison cf computed and recorded data 
for injector pistor displacement. As can be seen, the 
agreement is very good. This is not very significant in 
that all models tried, as well as the original analog 
computer model, were able to correctly predict injector 
pistcn metier. 

Figure (5) is a comparison of analytical results and 
experimental cata for breech chamber pressure. The computer 
model fellows the shape of the experimental curve fer the 
duration cf the time that the projectile is in motion (0-30 
msec) . It oscillates very rapidly but does not fall eff to 
a 2 erc value at the end of the projectile motion. The 
freguency of these initial oscillations is approximately 1.0 
khz which is within the -frequency range of the Kaman 
pressure transducers used in testing. However, the mounting 
of these transducers within a connecting cavity instead of 
flush with the chamber wall could have led to them 
experiencing a reduced, lagging freguency response. It is 
felt that seme, if not all, of the computed pressure 
oscillations must exist as evidenced by the close 
correspondence of the first two peaks in Fig. (5) . The 
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accuracy cf the experimental pressure reading could also be 
considered as being reduced by the location of the pressure 
transducer in the test setup. It is possible that the 
monitoring pressure tap was too close to the end cf the 
breech chamber and may not have sensed the full chamber 
pressure during the latter part of the projectile's travel 
when the fluid velocity is greatest. It is felt that to 
adeguately describe the pressure decrease toward zero, it 
would be necessary to completely account for exact changes 
in system geometry. To do this would require using a 
distributed parameter approach featuring variable lengths 
fcr the injector chamber and breech chamber. 

The pressure transients which occur after the projectile 
steps are ccnplex interference phenomena which are not fully 
described by the model. However, the peak pressure and the 
natural freguency of oscillation of the computer model 
compare favorably with the experimental data. 
Unfortunately, the damping characteristics cf the model do 
not follow the experimental data well. The system 
parameters which determine these values will be discussed in 
the next section. 

It should be ~ reiterated that the lumped parameter 
approach is at approximation,. In a distributed parameter 
approach, pressure and velocity would vary continuously with 
distance within the system as well as with time. In the 
lumped parameter approach, these changes are assumed to 
occur only at the input and output of a building block; 
therefore, pressure and velocity are considered constant 
within the building block. This will always lead to some 
discrepancies when comparing model data to experimental data 
taken at at a fix ed point. 

Figure (6) shows the simultaneous pressure history for 
several areas of the LPG feed system. The input ram 
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pressure , injector chamber pressure, connecting line 
pressure, ard breech chamber pressure are graphed as 
furcticns cf time. 

Figures (7) and (8) show the simultaneous velocities at 
several points as predicted by the computer model. The 
irjectcr piston velocity, fluid propellant velocities at the 
beginning aid end of the connecting line, and projectile 
velocity are plotted as functions of time. 

The injector and projectile displacements as functions 
of time are included as Figure (9) . As can be seen from 
Figures (7), (8), and (9), the movement of the injector 
piston is much mere rapid than that of the projectile. It 
is felt that this velocity difference is very much a 
furction cf system geometry and should not be considered as 
a generalization for all IPG feed systems. The NPS 
experimental apparatus being modeled had only small pressure 
drops between injector and breech chambers. In addition, 
the projectile slug was an order of magnitude lighter than 
the injectcr piston. In large scale systems, with 
significant pressure drops and very massive projectiles, it 
is guite possible that the injector pistcn‘s full stroke 
will occur significantly before the projectile has seated. 
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Figure 4 - COMPARISON OF INJECTOR PISTON DISPLACEMENT 



-1000 kPa 




(l-OL x ISd) 9anss9wi<j aaquieqo 



31 



Figure 5 - COMPARISON OF BREECH CHAMBER PRESSURE 
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Time (Sec. x 10^) 

Figure 6 - SIMULTANEOUS PRESSURE HISTORIES 
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Figure 7 - SIMULTANEOUS INJECTOR 
VELOCITY HISTORIES 




(^01 * oas/UL) /C1L3019A 



34 



Figure 8 - SIMULTANEOUS BREECH CHAMBER VELOCITY HISTORIES 
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VI 



ANALYSIS of system performance 



A. PARAKETEE S AFFECTING PROJECTILE RAM TIME 

Several parameters can be varied to "tune” the model to 
correct ram time. In equations (1) and (9) (the force 
balances cn the injector piston and projectile) empirical 
constants are included which account for any static friction 
and back pressure (PDI = injector back pressure, PDS = 
projectile lack pressure). In equations (4) and (6), 
resistance coefficients are introduced to account for system 
resistances. Finally in equations (2), (5), and (8), 
prcpellant effective bulk modulus values appear. 

The value for the injector and projectile total back 
pressure is difficult to determine precisely. Reference 1 
cited an experimentally determined value of 24 psi. The 
analog computer model presented in Ref. 2 used a value of 46 
psi for a ran pressure of 140 psi. The present model uses a 
value of 24 psi for PDI and 0.5 psi for PDS, or a total of 
24.5 psi tack pressure, which agrees with the experimental 
value. Ihe value of system back pressure will vary with 
driving pressure and the configuration of each LPG system's 
injector and breech chambers. 

Some tasic research has been conducted to relate the 
value of tie resistance coefficients to fluid properties. 
Onf ortunatelj these studies which are described in Refs. 4 
and 5, dealt with steady state fluid flow in constant 
diameter lines without flew restrictions , making their 
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results inapplicable to the present LPG feed system problem. 
The value cf EV^ can be related to the damping ratio of the 

"water hammer" pressure oscillations. The value of R?^ was 

determined bj an iterative process to obtain the best fit to 



experimental data. Figure 

ram-time on RV and RV , 
1 2 



(10) shows the dependence of 
keeping all ether variables 



constant. is expected intuitively, increasing the fluid 



resistance creates a larger system pressure drop, resulting 
in a lower pressure exerted on the face cf the projectile 
and hence slower ram times. The slope cf the two curves are 
almost identical sp that no significant advantage would be 
achieved by system designs which try to minimize either RV 



or RV to the detriment of the other. 
2 



It should be noted. 



however, that higher values for RV^ do tend tc slow down the 
ram time mere than high values for RV . Since the magnitude 



of E V ^ would probably depend mostly on the pressure drep at 
the gun valve, which seals the breech chamber, its design 



should he closely watched to ensure rapid ram times, 
likewise the system designer will have to pay close 
attention to the design of the piping and valves in a large 
scale LPG feed system to achieve optimum performance. 



The value cf the effective bulk modulus was determined 
from the uncamped natural freguency of the "water hammer" 
pressure cscillaticns as will be discussed in the next 
section. The bulk modulus is a fluid parameter which 
characterizes the spring effect of a liquid. The bulk 
modulus can be substantially lowered by the elasticity of 
the chamber and connecting line walls and the amount of 
entrapped gas present in the propellant fluid. As can be 



37 



seen iron Figure (11) , the value of the bulk modulus has 
only a small effect on ram time. Between bulk modulus 
values of 10,000 and 100,000 psi (the expected operating 
region of an operational system) the ram time is almost 
cccstant, varying less than a quarter of a millisecond. 
This is an encouraging result since minimizing entrapped gas 
is a difficult and cpstly design constraint. 
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Bulk Modulus (PSI 
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E. PARAMETEBS AFFECTING CHAMBER PRESSURE OSCILLATIONS 

As explained in Section III, after the projectile has 
stepped the chamber pressure can be described by a second 
order differential eguation (Egn. 10) . If the system 
described by this differential eguation is underdamped, the 
solution to the eguation is an exponentially damped sinusoid 
with a characteristic damping ratio and undamped natural 
f reguency . 

From Egn. (10), the system damping ratio, V is: 



If this daiiping ratio is increased, the peak chamber 
pressure as well as the dissipation period of the system 
pressure oscillations will be decreased. From a designer's 
point of ?iew, it would be advantageous for both of these 
guantities to be as low as possible. Hence, it would seem 
that an increase in system resistance to enhance damping 
would be desirable. However, as shown in the previous 
section, an increase in the system resistance also increases 
the ram time- an undesirable occurrence. Thus, alterations 
in system configuration which affect resistance coefficient 
values will have to be accomplished carefully to insure 
optimal sjstem performance. 

From Egn. (1C) # the undamped natural freguency cf the 
pressure oscillations in the LEG system is: 



The primary effect of ullage would be to decrease the value 




A3 RV2 

2\fpY, 
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of effective tulk modulus and consequently tc decrease the 
natural frequency of pressure oscillations. Hence, should 
it be found in future tests that the actual liquid 
propellant is highly sensitive to pressure oscillations, 
changes in ullage may be effected to decrease system 
ringing. Fortunately, as shown in the previous section, the 
effect of ullage variations on ram time is minimal. 
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VII. CONCLUSIONS 



Several important conclusions can be drawn from the 
output of the IPG ram feed computer model. Thus, it is 
realized that serious design consideration must be given to 
the values cf fluid resistances in the feed system in order 
tc minimize ram time and peak pressures, and to optimize the 
damping cf system pressure transients. Fortuantely, 
ccncommitant model results also show that the effect of 
entrapped gas is minimal for the overall operation of a LPG 
feed system. The model was unable to account fcr any 
cavitation effects so this is an area that must be a subject 
cf further studies. 

The fact that large pressure transients cccur during the 
prcjectile ramming cycle indicates the difficulties that 
might arise with proposed designs for liquid propellant guns 
with variable propellant volumes and projectile 
displacements. The pressure transients experienced in these 
systems will undoubtedly be very complex in nature. 

If further improvements are desired in the modeling of 
LPG feed systems, the next step should be the inclusion of 
thermody ramie effects. For any real fluid, resultant feed 
system pressure and velocity changes will cause temperature 
changes which may substantially affect the propellant fluid 
density. Consequently, propellant ignition characteristics 
at high rates cf fire will undoubtedly be affected. 

The lumped parameter LPG feed system computer model has 
been shown tc be in favorable agreement with experimental 
data. It is felt that this model has been sufficiently 
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validiated tc allow its use in more complicated LEG feed 
system designs. It is expected that this model will be used 
by the Naval Ordnance Station, Indian Head, Maryland to 
assist in studying the fluid dynamic characteristics of a 30 
mm scale model IPG feed system becoming operational in June 
1S76. It will be necessary to make the obvious changes to 
describe tie different geometry of the 30 mm scale model. 
The additional line lengths, and the larger pressure drops 
due to complex valve arrangements in the 30 mm system will 
reguire different values for model back pressure and 
resistance coefficients. The presence of an accumulator 
near the breech chamber will affect the value of the system 
camping ratic. It is possible that the different geometry 
and mass characteristics of the 30 mm model could result in 
the injector piston being slower than the projectile. 
However by incorporating all of these alterations, the 
manner in which each of the system parameters affect total 
system performance will be able to be predicted by the 
versatile lumped parameter computer model developed in this 
study. 

As such, computer simulation is a useful design tool. 
However, tc remain useful, it must be viewed in its proper 
perspective. Full scale prototype testing is the only 
conclusive method of demonstrating ordnance system 
performance. Unfortunately, testing is both costly and time 
consuming. Thus, computer simulation, no matter how simple 
or complex, can be used to designate critical testing 
instances and to identify those areas where design efforts 
will be most productive. In this way, the time involved in 
optimizing system performance as well as subseguent 
prcducticn ccsts can be reduced. In explosives research and 
ordnance design, which is still an empirical and necessarily 
hazardous science, computer simulation can be particularly 
useful in specifying the correct approach for design testing 
practices. In this respect, such models as the lumped 
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parameter model developed in this study can be considered as 
critical signposts at some of the many crossroads in systems 
development. However, it must be realized that as such they 
only point the way. The vehicle for arriving at an 
operational system can only be diligent research and 
engineering based on a measured progression of system 
demonstrations overseen by dedicated project managers aware 
of the many pitfalls along the way. 
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APPENDIX A 



COMPUTER PROGRAM LISTING 



The basis of the state variable method of systems 
analysis is the interpretation that the energy state of the 
components of a system completely describes the condition of 
a system. As defined in Ref. 6, the state of a system is the 
set of variables, the "state variables 1 ', which contain 
sufficient information about the present condition of the 
system tc permit the determination of all future time 
history cf the system - provided that all future inputs to 
the system are known. Therefore the energy state cf those 
elements which store energy completely describes the system. 
In the case of the lumped parameter model, these elements 
are the fluid inertia and capacitance and the inertia cf the 
injector pistcn and projectile. The energy state cf these 
storage elements can be described as a function of tine in 
terms of the state variables - pressure, velocity, and the 
system input pressure. Only those variables required to 
completely specify the state of the system need to be 
included . 

The first step in arriving at the computer program was 
to take eguaticns (1> through (9) and normalize, or solve 
them for the highest derivative. The state variables then 
become the pressure or velocity associated with these first 
derivatives. The computer program was developed by 
combining eguaticns and defining new variables as follows: 
X(1) = 1 
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x (2) = E 
X<3) = 1. 

X<4, = E * 

X<5) = V 
X (8) = T g 
X (9) = V g 

Equation (1) becomes: 

-Kfp Ai Ar Ai 

XD0T(1) X(l) - -- X(2) + -- Pr - -- Pd i 

Mp Mp Mp Mp 

Id the computer program/ PR is represented by X (3 0 ) and KFP 
and PDI become program constants/ C <7) and C (2 ) , 

respectively. The value of C(2) and C (7) are input on data 



cards. 



Using the fact that P^=P^ # eguation (2) becomes: 



XD0T(2) = p X(l) - - X(3) 

h L i 

where K is input to the computer program as C ( 3) . By 

combining eguaticns ,(3) and (4) and using the incompressible 

flew assumption A V =A V tire next eguation is derived: 

12 2 4 



XDOT (3) = 



Rho L 



( X(2) - Rv -i X(3) - X(4) ) 



where RV becomes C (4) . Next 
1 

incompressible flow assumptions 
yield: 



a 2 

eguation 



A V =A V 
12 2 4 



(5) 

and 



and 



the 



A V =A V 
2 5 3 6 



At K A, K 

XDOT (4) = — — X(3) - -2 - X(4) 

a 2 l 2 a 2 l 2 

and as before, C (3) =K. Eguations (6) and (7) combine with 

the flow velocity eguation A V =A V to yield: 

2 5 3 6 



XDOT (5) = 



A 3 



where RV 



_ ( x«, - Rv 2 j X(5) - X<8) ) 

becomes C (5) - Since V =V , eguation (8) becomes: 

6 7 



XDOT (8) = ~ X(5) - ~ X(9) 
L3 L3 
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where K again beccmes C (3 ) . Finally, equation (9) becomes: 

.(/pc A A 

XD0T(9) = X(9) + - X(8) - ^ Pds 

Ms Ms Ms 

where KFS and EDS become C (8) and C (9) respectively. The 
auxiliary eguaticns: 

XDOT (6) = X(l) ' 

XD0TC7) = X(9) 

are used tc calculate injector piston displacement, X (6) , 
and projectile displacement, X (7) . 
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